AFRL-RD-PS- 

TR-2014-0036 


AFRL-RD-PS- 

TR-2014-0036 


ADVANCED  CONTROL  FILTERING  AND  PREDICTION  FOR 
PHASED  ARRAYS  IN  DIRECTED  ENERGY  SYSTEMS 


James  Steve  Gibson 


University  of  California,  Los  Angeles 
Office  of  Research  Administration 
11000  Kinross  Avenue,  Ste  102 
Los  Angeles,  CA  90095-2000 


31  July  2014 


Final  Report 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  IS  UNLIMITED. 


AIR  FORCE  RESEARCH  LABORATORY 

Directed  Energy  Directorate 

3550  Aberdeen  Ave  SE 

AIR  FORCE  MATERIEL  COMMAND 

KIRTLAND  AIR  FORCE  BASE,  NM  87117-5776 


This  page  intentionally  left  blank. 


Approved  for  public  release;  distribution  unlimited. 


Notice  Page  for  Approved  for  Public  Release  TR  (Unlimited  Distribution) 


NOTICE  AND  SIGNATURE  PAGE 


Using  Government  drawings,  specifications,  or  other  data  included  in  this  document  for 
any  purpose  other  than  Government  procurement  does  not  in  any  way  obligate  the  U.S. 
Government.  The  fact  that  the  Government  formulated  or  supplied  the  drawings, 
specifications,  or  other  data  does  not  license  the  holder  or  any  other  person  or  corporation; 
or  convey  any  rights  or  permission  to  manufacture,  use,  or  sell  any  patented  invention  that 
may  relate  to  them. 

This  report  was  cleared  for  public  release  by  the  Air  Force  Research  Laboratory  [insert  TD 
site]  Public  Affairs  Office  and  is  available  to  the  general  public,  including  foreign  nationals. 
Copies  may  be  obtained  from  the  Defense  Technical  Information  Center  (DTIC) 
(http://www.dtic.mil). 


AFRL-RD-PS-TR-20 14-0036  HAS  BEEN  REVIEWED  AND  IS  APPROVED  FOR 
PUBLICATION  IN  ACCORDANCE  WITH  ASSIGNED  DISTRIBUTION 
STATEMENT. 


_ L /Signed/ /_ 

DAN  K.  MARKER,  DR-III 
Project  Officer 


//Signed// 

KENTON  T.  WOOD,  DR-IV,  DAF 
Chief,  Laser  Division 


This  report  is  published  in  the  interest  of  scientific  and  technical  information  exchange,  and  its 
publication  does  not  constitute  the  Government’s  approval  or  disapproval  of  its  ideas  or  findings. 

*Disseminated  copies  will  show  “//signature//”  stamped  or  typed  above  the  signature  blocks. 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the  data  needed, 
and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing  this  burden  to  Department  of 

Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that 
notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN 
YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

31-07-2014  Final  Renort 

3.  DATES  COVERED  (From  -  To) 

20121221  -  20140731 

4.  TITLE  AND  SUBTITLE 

ADVANCED  CONTROL  FILTERING  AND  PREDICTION  FOR  PHASED  ARRAYS  IN 
DIRECTED  ENERGY  SYSTEMS 

5a.  CONTRACT  NUMBER 

FA945 1-13-1-0265 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

62605F 

6.  AUTHOR(S) 

James  Steve  Gibson 

5d.  PROJECT  NUMBER 

5f.  WORK  UNIT  NUMBER 

D071 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

AND  ADDRESS(ES) 

University  of  California,  Los  Angeles 

Office  of  Research  Administration 

1100  Kinross  Avenue,  Ste  102 

Eos  Angeles.  CA  90095-2000 

8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory 

3550  Aberdeen  Ave  SE 

Kirtland  AFB,  NM  87117-5776 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

AFRL/RDLTS 

11.  SPONSOR/MONITOR’S  REPORT 

NUMBER(S) 

AFRL-RD-PS-TR-20 14-0036 

12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited.  377ABW-2014-1027;  Jan  8,  2015. 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

UCLA  research  under  this  project  has  developed  advanced  methods  for  control,  filtering,  prediction  and  system 
identification  in  adaptive  optics,  laser  beam  pointing  and  target  tracking.  These  methods  can  achieve  significant 
improvements  in  on-target  Strehl  ratios  and  tracking  jitter  for  phased-array  high  energy  laser  systems. 

15.  SUBJECT  TERMS 

Optics,  Strehl,  jitter,  phased-array,  high  energy  laser 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

SAR 

18.  NUMBER 
OF  PAGES 

17 

19a.  NAME  OF  RESPONSIBLE  PERSON 
Dan  K.  Marker 

a.  REPORT 

Unclassified 

b.  ABSTRACT 
Unclassified 

c.  THIS  PAGE 

Unclassified 

19b.  TELEPHONE  NUMBER  (include  area 
code) 

Standard  Form  298  (Rev.  8-98) 


1 


TABLE  OF  CONTENTS 


1.0  SUMMARY . 1 

2.0  INTRODUCTION . 2 

3.0  METHODS,  ASSUMPTIONS,  AND  PROCEDURES . 3 

3.1  Control  and  Filtering  Methods . 3 

3.2  Assumptions . 3 

3.3  Design,  Analysis  and  Simulation  Procedures . 3 

4.0  RESULTS  AND  DISCUSSION . 4 

4. 1  Feedback  Control  of  Phase  Ramp  with  Kalman  Prediction  from  Wrapped  Phase 

Measurements . 4 

4.2  Simulation  Model . 4 

4.3  Phase  Velocity  Range  for  Stable  Tracking . 5 

4.4  Effect  of  Sensor  Noise  and  Velocity  Variation  on  Closed-loop  Error  and  Velocity 

Estimates . 6 

5.0  CONCLUSIONS . 10 


Approved  for  public  release;  distribution  unlimited. 

iii 


LIST  OF  FIGURES 

Figure  1 :  SIMULINK  model  for  prediction  and  feedback  control  of  a  phase  ramp.  Mirror 


represented  by  integrator  with  sample  time  tsim.  The  model  shown  has  a  measurement 

delay  equal  to  the  length  of  two  control  sample  intervals . 5 

Figure  2:  Closed-loop  phase  error  12%.  Measurement  delay  =  2  control  sample  intervals. 

Velocity  vr  =  0.75  x27i/control  sample  interval.  Blue:  constant  velocity,  no  sensor 
noise.  Red:  velocity  standard  deviation  =  0.1  x27t,  sensor  noise  standard  deviation  = 

0. 1  x27i.  Phase  error  sampled  at  the  fast  simulation  rate . 6 

Figure  3:  Closed-loop  velocity  estimate  12%.  Measurement  delay  =  2  control  sample  intervals. 
Velocity  vr  =  0.75  x27i/control  sample  interval.  Blue:  constant  velocity,  no  sensor 
noise.  Red:  velocity  standard  deviation  =  0. 1  x27t,  sensor  noise  standard  deviation  = 

0. 1  x27i.  Phase  error  sampled  at  the  fast  simulation  rate . 7 

Figure  4:  Closed-loop  phase  error  /  2%.  Measurement  delay  =  2  control  sample  intervals. 

Velocity  vr  =  0.75  x27i/control  sample  interval.  Blue:  constant  velocity,  no  sensor 
noise.  Red:  velocity  standard  deviation  =  0.1  x2jt,  no  sensor  noise.  Phase  error  sampled 

at  the  fast  simulation  rate . 8 

Figure  5:  Closed-loop  velocity  estimate  /  2%.  Measurement  delay  -  2  control  sample  intervals. 
Velocity  vr  =  0.75  x27i/control  sample  interval.  Blue:  constant  velocity,  no  sensor 
noise.  Red:  velocity  standard  deviation  =  0.1  ^2%,  no  sensor  noise.  Phase  error  sampled 
at  the  fast  simulation  rate . 9 


Approved  for  public  release;  distribution  unlimited 

iv 


1.0  SUMMARY 


University  of  California,  Los  Angeles  (UCLA)  research  under  this  project  has  developed 
advanced  methods  for  control,  filtering,  prediction  and  system  identification  in  adaptive  optics, 
laser  beam  pointing  and  target  tracking.  These  methods  can  achieve  significant  improvements  in 
on-target  Strehl  ratios  and  tracking  jitter  for  phased-array  high  energy  laser  systems.  The  main 
technical  contribution  of  the  research  under  this  grant  is  a  method  for  phase  unwrapping, 
prediction  and  correction  in  phased  arrays.  The  new  method  combines  phase  unwrapping  logic 
with  Kalman  filtering  and  prediction  in  a  feedback  control  loop.  Analysis  and  simulation  have 
shown  that  the  new  method  can  track  and  correct  optical  phase  with  slowly  varying  phase 
velocity  and  that  the  method  is  reasonably  robust  to  sensor  noise.  Following  is  a  detailed 
discussion  of  this  method  and  representative  simulation  results. 
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2.0  INTRODUCTION 


Under  this  grant,  research  at  UCLA  has  developed  advanced  methods  for  control,  filtering, 
prediction  and  system  identification  in  adaptive  optics,  laser  beam  pointing  and  target  tracking. 
These  methods  can  achieve  significant  improvements  in  on- target  Strehl  ratios  and  tracking  jitter 
for  phased-array  high  energy  laser  systems. 

A  primary  objective  of  the  proposed  research  is  to  develop  novel  methods  for  wavefront 
prediction  and  correction  and  reduction  of  tilt  jitter,  such  as  two  methods  that  increase 
bandwidths  of  wavefront  correction  and  jitter  control.  An  equally  important  objective  will  be  to 
develop  analysis  and  simulation  tools  to  evaluate  closed-loop  performance  of  phased  array 
systems  with  a  variety  of  control  system  designs  and  subject  to  realistic  disturbances, 
measurement  and  reconstruction  errors  and  system  latencies. 

The  main  technical  contribution  of  the  research  under  this  grant  is  a  method  for  phase 
unwrapping,  prediction  and  correction  in  phased  arrays.  The  new  method  combines  phase 
unwrapping  logic  with  Kalman  filtering  and  prediction  in  a  feedback  control  loop.  Analysis  and 
simulation  have  shown  that  the  new  method  can  track  and  correct  optical  phase  with  slowly 
varying  phase  velocity  and  that  the  method  is  reasonably  robust  to  sensor  noise.  Following  is  a 
detailed  discussion  of  this  method  and  representative  simulation  results. 
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3.0  METHODS,  ASSUMPTIONS,  AND  PROCEDURES 

3.1  Control  and  Filtering  Methods 

The  primary  methods  used  to  design  the  predictor  and  feedback  control  loop  are  the 
following: 

•  Kalman  filtering  and  prediction 

•  Feedback  control 

•  Phase  Unwrapping 

3.2  Assumptions 

It  is  assumed  that  a  possibly  noisy  measurement  of  wrapped  optical  phase  is  available  to 
the  digital  feedback  control  system. 

3.3  Design,  Analysis  and  Simulation  Procedures 

The  research  procedure  consists  of  designing  the  predictor,  phase  unwrapping  logic  and 
feedback  control  law,  analyzing  the  closed-loop  system  theoretically  to  guarantee  stability  and 
simulating  the  closed-loop  system  in  SIMULINK. 
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4.0  RESULTS  AND  DISCUSSION 

4.1  Feedback  Control  of  Phase  Ramp  with  Kalman  Prediction  from  Wrapped  Phase 
Measurements 

This  discussion  summarizes  results  from  a  study  of  mod  47t  phase  measurement  and  the 
effects  of  sensor  noise  and  velocity  variation  on  the  estimation,  prediction  and  control  of  a  phase 
ramp.  Figure  1  shows  a  closed-loop  system  with  two  rates,  one  rate  for  simulation  and  the 
second  rate  for  the  controller.  The  mirror  is  represented  by  an  integrator,  so  that  the  control 
command  is  a  rate  command  to  the  mirror.  This  model  assumes  that  the  position  of  the  mirror 
can  be  measured  and  used  by  the  control  loop,  although  this  measurement  probably  is  not 
necessary. 

4.2  Simulation  Model 

There  are  three  differences  between  the  current  SIMULINK  model  and  the  previous 
models:  The  measurement  is  mod  4ji  in  the  current  model  instead  of  mod  2n  in  the  previous 
models,  white  sensor  noise  is  added  to  the  phase  measurement  before  the  47t  block,  and  the 
current  model  provides  for  a  band-limited  disturbance  to  be  added  to  the  phase  velocity.  The 
results  discussed  here  were  generated  with  only  the  white  sensor  noise.  The  study  with  velocity 
disturbance  is  in  progress  and  will  be  reported  later. 

The  controller  consists  of  the  blue  blocks  in  Figure  1.  The  open-loop  phase  is  the  signal 
to  be  predicted  and  corrected.  The  Kalman  predictor  and  unwrapping  function  take  the  delayed 
closed-loop  phase  error  (mod  47t)  as  the  only  input,  and  generate  a  prediction  of  the  open-loop 
phase.  This  prediction  is  subtracted  from  the  true  open-loop  phase  to  form  the  closed-loop  phase 
error. 

The  sample-and-hold  rate  for  control  is  defined  to  be  tcontr  =  1 .  The  simulation  rate, 
denoted  by  tsim,  normally  is  some  fraction  of  tcontr.  For  the  simulation  results  in  the  subsequent 
figures,  tsim  =  1/40,  and  the  measurement  delay  is  either  one  control  sample  interval  (=  40  tsim ) 
or  two  control  samples  (=  40  tsim).  A  few  simulations  with  tsim  =  1/200  yielded  essentially 
identical  results  to  those  presented  here,  so  tsim  =  1/40  was  used  for  most  of  the  numerous 
simulations  of  which  the  results  here  are  representative. 
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4.3  Phase  Velocity  Range  for  Stable  Tracking 

The  most  important  parameter  is  the  ratio  of  the  phase  velocity  vr  to  the  control  sample 
interval.  Since  that  sample  interval  is  1  in  this  model,  the  important  ratio  is  vr.  For  the  mod  4n 
phase  measurement,  the  current  controller  can  track  the  phase  and  produce  a  reliably  stable 
closed-loop  error  for  phase  velocities  with  magnitude  up  to  0.75><27i  /  control  sample  period. 
The  closed-loop  track  loop  is  stable  for  somewhat  larger  velocities  for  sufficiently  small  sensor 
noise  and  some  values  of  the  initial  phase.  This  suggests  that  further  refinement  of  the  Kalman 
predictor  and  control  loop  might  produce  a  stable  track  loop  for  velocities  somewhat  larger  than 
0.75x271,  but  possibly  not. 


Multi-Rate  Discrete  Time  Simulation 


Discrete 


Figure  1.  SIMULINK  model  for  prediction  and  feedback  control  of  a  phase  ramp.  Mirror 
represented  by  integrator  with  sample  time  tsim.  The  model  shown  has  a  measurement  delay 
equal  to  the  length  of  two  control  sample  intervals. 
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4.4  Effect  of  Sensor  Noise  and  Velocity  Variation  on  Closed-loop  Error  and  Velocity 
Estimates 

Figures  2  and  3  show,  respectively,  the  closed-loop  phase  error  and  velocity  estimate  for 
vr  =  0.75x271  /  control  sample  period  (approximately  the  maximum  velocity  that  this  controller 
can  track  with  reliable  stability)  and  zero-mean  white  sensor  noise.  The  measurement  delay  was 
two  control  sample  periods  and  the  sensor  noise  and  velocity  both  had  standard  deviation  = 

0. 1  x27i.  (The  standard  deviation  of  the  sensor  noise  equals  the  RMS  value,  since  the  noise  has 
zero  mean.)  For  Figures  4  and  5,  the  velocity  had  standard  deviation  =  0. 1  x27t  but  there  was  no 
sensor  noise. 


Closed-loop  Phase  Error/2n 


Figure  2.  Closed-loop  phase  error  /  27t.  Measurement  delay  =  2  control  sample  intervals. 
Velocity  vr  =  0.75x27i  /  control  sample  interval.  Blue:  constant  velocity,  no  sensor  noise.  Red: 
velocity  standard  deviation  =  0.1  x27t,  sensor  noise  standard  deviation  =  0.1  x27i.  Phase  error 

sampled  at  the  fast  simulation  rate. 
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Closed-loop  Velocity  Estimate/271 


Figure  3.  Closed-loop  velocity  estimate  /  2n.  Measurement  delay  =  2  control  sample  intervals. 
Velocity  vr  =  0.75x27t  /  control  sample  interval.  Blue:  constant  velocity,  no  sensor  noise.  Red: 
velocity  standard  deviation  =  0.1  x2ti,  sensor  noise  standard  deviation  =  0.1  x2tt.  Phase  error 

sampled  at  the  fast  simulation  rate. 
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Closed-loop  Phase  Error/2n 


Figure  4.  Closed-loop  phase  error  /  2n.  Measurement  delay  =  2  control  sample  intervals. 
Velocity  vr  =  0.75x2jt  /  control  sample  interval.  Blue:  constant  velocity,  no  sensor  noise.  Red: 
velocity  standard  deviation  =  0.1x27i,  no  sensor  noise.  Phase  error  sampled  at  the  fast  simulation 

rate. 
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Closed-loop  Velocity  Estimate/271 


Figure  5.  Closed-loop  velocity  estimate  /  2k.  Measurement  delay  =  2  control  sample  intervals. 
Velocity  vr  =  0.75x2jt  /  control  sample  interval.  Blue:  constant  velocity,  no  sensor  noise.  Red: 
velocity  standard  deviation  =  0. 1  x2tt,  no  sensor  noise.  Phase  error  sampled  at  the  fast  simulation 

rate. 
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5.0  CONCLUSIONS 

This  project  has  developed  advanced  methods  for  control,  fdtering,  prediction  and  system 
identification  in  adaptive  optics,  laser  beam  pointing  and  target  tracking.  These  methods  can 
achieve  significant  improvements  in  on-target  Strehl  ratios  and  tracking  jitter  for  phased-array 
high  energy  laser  systems.  The  main  technical  contribution  of  the  research  under  this  grant  is  a 
method  for  phase  unwrapping,  prediction  and  correction  in  phased  arrays.  The  new  method 
combines  phase  unwrapping  logic  with  Kalman  filtering  and  prediction  in  a  feedback  control 
loop.  Analysis  and  simulation  have  shown  that  the  new  method  can  track  and  correct  optical 
phase  with  slowly  varying  phase  velocity  and  that  the  method  is  reasonably  robust  to  sensor 
noise.  Further  research  should  extend  the  methods  developed  here  to  multiple  apertures  and 
higher-order  wavefront  correction  in  phased-array  systems. 
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